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A new formulation of the equation for calculation of air density has been developed. The Cohen and Taylor 
value of the gas constant, currently accepted values of the atomic weights, and recent determinations of 
abundances of the various constituents of air have been used. The abundance of carbon dioxide has been treated 
as a variable and a factor enabling convenienl adjustment of the apparent molecular weight of air fpr deviation of 
carbon dioxide abundance from a background value has been derived. A new table of the compressibility factor for 
the range of pressure and temperature of interest in standards laboratories has been calculated using reeeiitU 

determined values of virial coefficients. The enhancement factor, which has usually been ignored, lias been 

explicitly included. A simple equation for the calculation of enhancement factor has been lilted to values in the 
range of |>re>snre and temperature of interest. A simple equation for the calculation of saturation water vapor 
pressure has been fitted. Uncertainties, random and systematic, in the parameters and in the measurement of 
environmental variables and consequent uncertainties in calculated ait density have been estimated. 

Application of the equation to air buoyancy determination and the transfer of the mass unit at the various 
national standards laboratories has been made. 
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1. Introduction 

The transfer of the mass value from one object, such as 
the International Prototype Kilogram, to another object is 
accomplished 1)\ comparison of the objects bv means of a 
balance. The difference in buoyant force on the two objects 
is proportional to the difference in their displacement vol- 
umes and to the air density. The air density is conventionally 
calculated using; an equation based on the equatioti of state 
of an air-water vapor mixture. A new formulation of the air 
density equation is developed below. 

2. Development of the Air Density Equation 

The total pressure, P, the total volume, V, and the 
absolute temperature, 7\ of a mixture of ideal gases are 
related by the ideal gas equation, 



PV = nR7\ 



(1) 



where n is the number of moles of the mixture and R is the 
universal gas constant. In terms of density, p, rather than 
volume, eq (1) becomes 



P = P ^ (2, 



vhere M is the apparent molecular weight of the mixture. 



For a mixture of dry air (indicated bv subscript a ) and 
water vapor (subscript w), p and VI are the density and appar- 
ent molecular weighl respectivel) of the air-water vapor mix- 
ture. Since 



M 



m m (l -f- m„ 



n n a + n w 



(3) 



where m is the mass of the mixture and // is the number of 
moles of the mixture, 




M = n a M a + n u M w =M 



By introducing the water vapor mixing ratio, r: 
mass of water vapor n w M w 



(4) 



mass of dry air n a M a " 



M , 
and by designating the ratio by e, eq (4) becomes 

lVl n 



(5) 



M = 



Mail + r) 



1 



(6) 
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Substituting (6) in (2) and noting that the effective vapor 
pressure, e', of water in moist air is defined [l] 1 by: 



(e + r) 



P, 



then 



R 
P = p — T 

M„ 



1 + (e - 1) 



(7) 



(8) 



Equation (8) is the ideal gas equation for a mixture of dry air 
and water vapor with a vapor pressure of e' . If the air-water 
vapor mixture behaved as a mixture of ideal gases, 
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where Z is the compressibility factor. Since the mixture is 
not ideal the magnitude of the non-ideality is reflected in the 
departure of Z from 1 and (9) becomes 



R 

P = p — TZ 
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1 + (€ - 1) 



PA 



(10) 



Equation (10) is the real gas equation for a mixture of dry air 
and water vapor. By rearrangement of eq (10), the expression 
for the air density is 
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3. Specification of the Values of the 
Parameters in the Air Density Eq (1 1) 

3.1 . Universal Gas Constant, R 

The value of the molar gas constant, /?, listed in a 
compilation by Cohen and Taylor [2], is 8.31441 ± 0.00026 
JK _1 mol -1 . Recently, Quinn et al. [3] made a new determi- 
nation of R by measuring the speed of sound in argon by 
means of an acoustic interferometer. Their value was 
8315.73 ± 0. 17 JK" 1 kmol" 1 . Gammon [4] recently deduced 
a value of R from measurements of the speed of sound in 
helium; his latest reported value is 8315.31 ± 0.35 JK -1 
kmol -1 [5], which is in close agreement with the Quinn et al. 
value. Rowlinson and Tildesley [6] have recently interpreted 



the experimental measurements of Quinn et al. and arrived 
at a value of 8314.8 ± 0.3 JK -1 kmol" 1 , which is in 
agreement with the Cohen and Taylor value within the 
uncertainties assigned to the values. 

We choose at present to use the Cohen and Taylor value 
with the realization that in the future it might be replaced by 
a new value. 



3.2. Apparent Molecular Weight of Air, M a 



The apparent molecular weight of dry air, M a , is calcu- 
lated using the relationship 



M a = 2 M iXi 



(12) 



1 Figures in brackets indicate references at the end of this paper. 



where each M t is the molecular weight of an individual 
constituent and x t is the corresponding mole fraction. The 
molecular weights and typical mole fractions of the constitu- 
ents of dry air are tabulated in table 1. Other constituents 
are present in abundances which are negligible for the 
present application. 

The values of the atomic weights of the elements are taken 
from reference [7] and are based on the carbon-12 scale. The 
molecular weights are taken to be the sums of the atomic 
weights of the appropriate elements. 

The value for the abundance of oxygen is taken from 
reference [8]. The value for the abundance of carbon dioxide 
is taken from a recent unpublished compilation of data on 
atmospheric concentration of carbon dioxide at seven loca- 
tions throughout the world. It must be emphasized that 
0.00033 is the mole fraction of CO2 in the atmosphere and 
should be considered to be a "background" value. The mole 
fraction of CO2 in laboratories, which is of course the value 
of interest here, is in general greater than 0.00033 and is 
variable. For example, three samples of air taken from a 
glove box in the Mass Laboratory at NBS had a mean value of 
0.00043, and four samples of laboratory air taken at the 
National Center for Atmospheric Research in Boulder, Colo- 
rado had a mean value of 0.00080. Clearly, then, the 
optimum utilization of the air density calculation would 
necessitate a measurement of CO2 abundance on an air 
sample taken at the time of the mass comparison. 

One of the options one has in dealing with the variability 
of CO2 abundance is to select a reference level (for example, 
0.00033 or 0.00043) and to provide an adjustment to M a to 
account for known departures from the reference level. 
Gluekauf [9], in discussing the variation of the abundance of 
oxygen in the atmosphere, stated that "all major variations of 
the O2 content must result from the combustion of fuel, from 
the respiratory exchange of organisms, or from the assimila- 
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tion of CO2 in plants. The first process does not result in 
more than local changes of O2 content, while the latter two 
processes, though locally altering the CO2/O2 ratio, leave 
their sum unchanged." The assumed constancy of the sum of 
the 2 and CO2 abundance simplifies the adjustment of M a 
to account for departures from the CO2 reference level and 
simplifies the estimation of the uncertainty in air density due 
to an uncertainty in CO2 abundance. The constancy of the 
sum is expressed by the equation (for convenience, the 
subscript i has been replaced by the chemical symbol): 



*co 2 "I" Xq 2 - cons 



nstant = 0.20979. 



(13) 



The contribution of O2 and C() 2 to the apparent molecular 
weight of dry air is 

M 02 x 0i + M COt x COs = 31.9988 x 0i + 44.0098 x C02 . ( 14) 



From (13), 



and 



s = 0.20979 - *co a 



M x ()2 + M C o 2 3 C o 2 
Therefore, 



12.01 



• v co.. 



(15) 



+ 0.7130. (16) 



volume concentration in U.S. Standard Atmosphere, 1976 

[in. 

From the data of table 1, the apparent molecular weight of 
dry air with a CO2 mole fraction of 0.00033 is calculated by 
(12) to be 28.963. 

TABLK /. Composition of <h\ on 



Constituent 


Abundance 
(mole fraction) 


Molecular 
Weight 


Nitrogen (N 2 ) 


0.78102 


28.0134 


Oxygen (0 2 ) 


.20946 


31.9988 


Carbon Dioxide (C0 2 ) 


.00033 


44.0098 


Argon (A) 


.00916 


39.948 


Neon (Ne) 


.00001818 


20.179 


Helium (He) 


.00000524 


4.00200 


Krypton (Kr) 


.000001 11 


83.80 


Xenon i \<i 


.000000087 


131.30 


Hydrogen (H 2 ) 


.0000005 


2.0138 


Methane (CH4) 


.000001.") 


10.0426 


Nitrous Oxide (N 2 0) 


.0000003 


44.0128 



3.3. Compressibility Factor, Z 

The compressibility factor is computed using the virial 
equation of slate of an air-water vapor mixture expressed as 
a power series in reciprocal molar volume. 



8(M a ) = 8 [M ,*o, + M co ,*co 2 ] 
= 12.011 8(* co ,), 



(17) 



that is, the variation in M a due to a variation in (X) 2 
abundance is equal to 12.01 1 (the atomic weight of carbon) 
multiplied by the variation in C0 2 abundance. The variation 
in M a due to the difference between the reference levels 
0.00033 and 0.00043 is thus 0.0012 g mol" 1 which corre- 
sponds to a relative variation of 41 p. p.m. in M a and a 
corresponding relative variation of 41 p.p.m. in air density. 

The adjusted M a accounting for the departure of the C0 2 
abundance from the reference level of 0.00033 becomes 



M fl 



M afm + 12.011 |> C o 2 - 0.00033], (18) 



where M a is the apparent molecular weight of dry air with 
a CO2 mole fraction of 0.00033. 

The value of the abundance of argon in dry air, 0.00916, 
is that calculated from the mass spectrometric determination 
of the ratio of argon to argon and nitrogen b) Hughes [ 1()|. 

The value for the abundance of nitrogen was arrived at by 
the usual practice of inferring nitrogen abundance to be the 
difference between unity and the sum of the mole fractions of 
the other constituents. 

The abundances of the constituents neon through nitrous 
oxide in table 1 were taken to be equal to the parts per 



Z 



El 

RT 



+ 



Bmix , Ci 



(19) 



and expressed as a power series m pressure, 



Pv 



Z =W= ] + / ''mix/' + ^mi*/' 2 + (20) 

A/ 



the molar 



ne, H v 



and II' m i X are second virial 



where v is trie molar volum* , u m [ X 
coefficients and C mix and C' mix are third virial coefficients 
for the mixture. The virial coefficients of the pressure series 
are related to the virial coefficients of the volume power 
series by 



D mix 



RT 



mix 



(RTf 



(21) 



(22) 



Each mixture virial coefficient is a function of the mole 
fractions of the individual constituents and the virial coeffi- 
cients for the constituents. The latter are functions of 
temperature only. 
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The second interaction (cross) virial coefficient of moist 
air, B aw , is one of the contributors to B' mix and expresses 
the effects of interaction between an air molecule and a water 
molecule. The values oiB aw used in the calculation of Z are 
experimental values which strictly apply to C02-free air. 

Using (28) and (29), below, and the virial coefficients [12, 
16J provided by Hyland,, a table of compressibility factor, 
Z, for C02-free air, table 2, has been generated. 

In the absence of values for virial coeffieients for air-C0 2 
mixtures in the temperature range of interest, the effect of 
the variability of C0 2 abundance on Z has been estimated to 
be of the order of 1 ppm and, therefore, negligible. Also, the 
effect of the variability of C0 2 abundance on Z due to the 
interaction of C0 2 with water vapor has been estimated to be 
of the order of 0.1 ppm and, therefore, negligible. Conse- 
quently, table 2 is applicable to moist air containing reason- 
able amounts of C0 2 . 

For temperatures and/or pressures outside the range of 
table 2, the table of compressibility factor of moist air (also 
C02-free) in the Smithsonian Meteorological Tables [15] can 
be used, with some loss of precision since the listing there is 
to the fourth decimal place. 

3.4. Ratio of the Molecular Weight of Water to the 
Molecular Weight of Dry Air 

The molecular weight of water is 18.0152 [7]. The ratio, 
e, of the molecular weight of water to that of dry air is, 
therefore, 0.62201 for dry air with a C0 2 mole fraction of 
0.00033. 

4. Uncertainty in Calculation of Air Density 

In estimating uncertainties we shall report them as 1 
standard deviation and we shall follow the suggested practice 
of Eisenhart [13, 14] in stating separately the random and 
systematic components. 

4.1. Uncertainty in R 

The uncertainty in the value of the molar gas constant is 
that reported by Cohen and Taylor [2], ± 0.00026 JKT 1 
mol -1 . The corresponding random relative uncertainty in p is 
3.1 X 10" 5 . 

4.2. Uncertainties in M a 

The uncertainties in M a will be taken from the estimated 
uncertainties in particular measurements of the abundances 
of the constituents of air. The uncertainty in the 2 abun- 
dance (mole fraction) of air [8] is separable into a random 
component and a systematic component. The random uncer- 
tainty is estimated to be ±0.00001, the systematic compo- 
nent is ±0.00006. 
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The random component of the uncertainty in the argon 
abundance is inferred from the precision of Hughes 1 mea- 
surements [10] to be ±0.00003. No assignmenl of systematic 
uncertainty was made by Hughes. 

The uncertainty in the CO2 abundance has been mentioned 
earlier with respect to the variability of the CO2 abundance 
in the laboratory. For a sample of air taken at the time a 
mass comparison is made, the estimated uncertainty in the 
subsequent mass spectrometric determination of CO2 abun- 
dance is ±0.00003 at the 0.00033 level. Since the measure- 
ments made by the mass spectrometric method are consid- 
ered to be very precise, the estimated uncertainty is consid- 
ered to be systematic. 

Since the N 2 abundance is the difference between unity 
and the sum of the mole fractions of the other constituents. 
the random component of the uncertainty in the N 2 abun- 
dance is found l)\ combining by quadrature the random 
components of the uncertainties of the other three major 
components to be ± 1 X 10 \ For the systematic component, 
however, the practice ol finding the V> abundance b\ 
difference provides for very significant reduction of uncer- 
tainty. This is shown in the following treatment. 

Equation ( 12) can be written: 

M (l = M 02 xq 2 + M C o 2 *co 2 

a- (23) 

+ M A x A + M n ./n 2 + 2 M> x i- 

1=5 

The mole fraction of N2 is, b\ difference, 

k 

x N = 1 - x - .in, - x A - 2j x i- (24) 

1=5 

If we ignore the last term in (28), differentiate and go to 
finite differences, 



MM (l ) Xi = M 02 Axo 2 + Af C o 2 A*co 2 
+ M A Ax A + M N2 
(-A*q 2 - Aa C02 - Ax A ) 
= (M 0i - M N2 )A^o 2 

+ (M C o 2 - Mn 2 )A.v ( , )2 
+ (M A - M N jA.v.,. 



found to be ±4 X L0 by combining by quadrature the 
product of the molecular weight and the random component 
of uncertainty in abundance for each of the four major 
constituents and dividing l>\ \l a . The uncertainty in M a due 
to the uncertaint) in the value of the atomic or molecular 
weights of the /-th constituent is 



A(Af„) A/| . = Xi&Mt 



(26) 



(25) 



With the substitution of the appropriate systematic uncertain- 
ties and the molecular weights into (25), the systematic 
components of the relative uncertainty in M Q due to uncer- 
tainties in abundance of the constituents if found to be ±3 
X 10 '\ The random component of the relative uncertainty is 



The uncertainties in the values of the atomic or molecular 
weights are inferred from reference [7|.~ For the four major 
constituents they are: for 2 , ±0.0002; for C0 2 , ±0.0005; 
for ,4, ± 0.001; and for N 2 , ±0.0001. These uncertainties 
are considered to be systematic. The sum of the four terms 
represented l>\ (26) is ±0.0001. corresponding to a relative 
uncertaint) in M a of ±4 x 10 ,; . 

The random component of the overall relative uncertainty 
in M a is, by quadrature, ±4 x 10 \ The systematic 
component ol the overall relative uncertaint) in M a is ±3 x 
10 , the sum of the component due to the systematic 
uncertainty in abundance and that due to the uncertainty in 
atomic or molecular weight. The corresponding random and 
systematic components of the relative uncertainty in f) are 
±4 X 10 " 5 and ±3 X 10 \ respectively. 

4.3. Uncertainty in Z 

The uncertaint) in Z, the compressibilit) factor, is esti- 
mated from the various uncertainties in Z due to the 
uncertainties in the virial coefficients | 12, 16]. The estimated 
relative uncertainty in/ at 293.15K, 10132.") Pa (1 atmos- 
phere) and SO percent relative humidit) is ±1.7 x 10 -,> . 
The corresponding relative uncertainty, taken to be sys- 
tematic, in p is ±1.7 X L0~ 5 . 

4.4. Uncertainty in e 

The uncertainty in M w , the molecular weight of water, is 
±0.0005 [7] and is treated here as systematic. The random 
component of the uncertainty in e, the ratio ofM w to M a , is 
calculated to be ±1 X 10 -5 and the systematic component is 
estimated to be ±8 X 10~ 6 . The corresponding uncertainties 
in the factor in (11) involving 6, 



1 + (6 



(27) 



at 293. L5K, 101325 Pa and 50 percent relative humidity are 
±1 X 10~ 7 and ±9 X 10~ 8 , respectively. The corresponding 
relative uncertainties in the term and consequently in p are 
±1 X 10" 7 and ±9 X 10~ 8 , respectively. 



2 The uncertainty assigned in the reference is divided by 3 to estimate the 
uncertainty at 1 standard deviation. 
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4.5. Combined Relative Uncertainties in p Due to the 
Uncertainties in M a , R, Z and e 

The random and systematic components of the relative 
uncertainty in p at 293. 15K, 101325 Pa and 50 percent 
relative humidity, due to the uncertainties in A/ a , R, Z and 
e are tabulated in table 3. 



TABLE 3. Components of the relative uncertainty in p due to uncertainties 
M a , R, Z and e 



Parameter 


Random 


Systematic 


Ma 


± 4 X 10" 5 


± 3 X 10~ 5 


R 
Z 

€ 


± 3.1 X JO" 5 




± 1.7 x 10" 5 
± 9 X 1(T 8 


± i x nr 7 


Combined 


± 5 x icr 5 


± 5 X 1(T 5 



The random components were combined by quadrature, 
the systematic components were combined by addition. 
These components of uncertainty can be considered to 
represent the uncertainty "■intrinsic" to the air density 
equation, that is, that which is contributed by the limitations 
on the present knowledge of M a , M w , R and Z. With these 
uncertainties must, of course, be combined the uncertainties 
in the environmental variables: P, T and relative humidity, 
and in the knowledge of the carbon dioxide abundance. 

5. Measurement of P, T and Relative 
Humidity 

In order to estimate the uncertainties in p due to uncer- 
tainties in the measurements of the environmental variables, 
P, T and relative humidity, we shall estimate the uncertain- 
ties in these measurements when made using the best 
applicable instrumentation and procedures. Therefore, the 
estimated uncertainties in p will be those contributed by the 
best possible measurements. 

5.1 . Pressure Measurement 

The state-of-the-art in pressure measurement [17] permits 
the measurement of the pressure in a laboratory with a 
random relative uncertainty of less than ±2 X 10 -4 . Calibra- 
tion of pressure measuring instruments against a primary 
standard of pressure contributes a systematic relative uncer- 
tainty of about ±3 X 10~ 5 . The corresponding Ap/p's 
contributed by uncertainties in the measurement of pressure 
in a laboratory in the vicinity of a balance case are ±2 X 
L0 _4 and ±3 x LO -5 . 

5.2. Temperature Measurement 

The measurement of temperature in a balance case, that is 
the temperature which determines the buoyant forces, is 



potentially the most critical measurement in terms of its 
effect on the uncertainty in the calculated air density. In the 
absence of experimental results, it is possible at this time to 
make only a rough estimate of the temperature uncertainty to 
be expected. If the balance case were instrumented with a 
network of thermopile junctions, for example, the measure- 
ments would be expected to have a standard deviation of 
about ±0.05K [18]. The systematic uncertainty is estimated 
to be of the order of ±0.01K. At a temperature of 293.15K, 
these uncertainties correspond to relative uncertainties of 
±2 X 10 -4 , and ±3 X 10~°, respectively. The correspond- 
ing Ap/p's are ±2 X 10~ 4 and ±3 X 10 -5 , respectively. 

5.3. Humidity Measurements 

The state-of-the-art in humidity measurement [19] permits 
the measurement of humidity in a balance case with a 
random uncertainty of =t0.5 percent relative humidity and a 
systematic uncertainty of ±0.3 percent relative humidity. 
These uncertainties correspond to relative uncertainties in 
the water vapor pressure factor |1 + (e — l)e'IP], in (11) 
and, therefore, to Ap/p, of ±4 X 10" 5 and ±3 X 10" 5 at 
293. 15K, 101325 Pa and 50 percent relative humidity. 

Since e' is the effective vapor pressure of water in moist 
mr, a word of caution with regard to inferring e from 
measurements of relative humidity is in order. Relative 
humidity, U, can be defined [20] by 



U = — X 100 percent, 



(28) 



where e s is the effective saturation vapor pressure of water in 
moist air. e s is greater than e s , the saturation vapor pressure 
of pure phase (i.e., water vapor without the admixture of air 
or any other substance) over a plane surface of pure ordinary 
liquid water, since the introduction of a second gas (air in 
this case) over the surface of the water increases the 
saturation concentration of water vapor above the surface of 
the water. This "enhancement" of water vapor pressure is 
expressed by the enhancement factor, /, which is defined by 



/= 



(29) 



The most recently published [16] experimentally derived 
value of/ at 293. 15K and 100000 Pa is 1.00400. Therefore, 
the common practice of inferring e' from measured U and 
tabulated values of e 8 introduces a significant error in e iff 
has been ignored. The corresponding relative error in p at 
293. 15K, 101325 Pa and 50 percent relative humidity is 
about 1.7 X 10" 5 . 

f is a function of temperature and pressure. In the present 
work, Hyland's values of/ [16] have been fitted to a three- 
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parameter equation in the pressure and temperature (/, °C) percent relative humidity, an error of 0.1 percent in calcu- 
ranges of interest in national standards laboratories. The lated air density results from an error of 0.29K in temperature 
resulting equation is measurement or a 101 Pa error in pressure measurement or 

an 11.3 percent relative humidity error in the measurement 
/= 1.00070 + 3.113 X \(r H P + 5.4 X IO" 7 * 2 . (30) of relative humidity. 



The expression fore' is found by combining (28) and (29) 
to be 



U 

e — fe<, 

!()()• ' 



(31) 



The systematic relative uncertainties in p due to the 
uncertainties assigned toy [16] and e s [22] are approximately 
± 1 X 10~ 6 and ±2 X \Q-\ respectively. 

For the temperature range of interest in the present 
application, any of several tables of e 8 , for example, refer- 
ences [21], [22], and [23], can be used. Besley and Bottom- 
ley [24] have recently published experimental values of e, in 
the temperature range 272.60 to 208. 04K. 

In the present work, the data of Besley and Bottomley in 
the temperature range 288.15 to 208. 04K and calculated 
values |22] for the remainder of the temperature range to 
301. 15K, have been fitted to a two-parameter equation. The 
resulting equation is 



e 8 = 1.7526 X 10 11 exp (-531 5.56/7'). 



(32) 



Values calculated using (32) are sufficiently (lose to 
experimental and calculated values, within ±0. 1 percent, to 
be used in the present application. 

5.4. Random and Systematic Components of the 

Relative Uncertainty in p Due to Anticipated 

Uncertainties in State-of-the-Art Measurements of P, T 

and Relative Humidity 

The random and systematic components of the relative 
uncertainty in p at 293. 15K, 101325 Pa and 50 percent 
relative humidity due to anticipated uncertainties in state-of- 
the-art measurements of the environmental variables are 
tabulated in table 4. The random components were combined 
by quadrature, the systematic components were combined by 
addition. 

It is clear from inspection of tables 3 and 4 that the 
dominant uncertainty in the calculation of air density is that 
contributed by the measurement of pressure, temperature 
and relative humidity, even when the measurements are 
made using the best instrumentation and procedures. Consid- 
erable care must therefore be taken in making measurements 
of these environmental variables in order to attempt to 
approach the precision and accuracy represented in table 4. 

To illustrate the effect of errors in the measurement of the 
environmental variables, at 293. 15K, 101325 Pa and 50 



Table 4. Components of the relative uncertainty in pdue to (Anticipated 
uncertainties in state-of-the-art measurements of P. T and relative humidity 



\ ariable 


Random 


Systematic 


P 

T 

Relative Humidit) 

Combined 


± 2 X 10 « 
±2X|()' 
± 4 x LO" 5 
± 3 x io 4 


± 3 x LO -5 
± 3 x L0" 5 

± 3 x io : > 
± 9 X 10 •> 



5.5. Carbon Dioxide Abundance 

As stated in section 3.2, the CA) 2 abundance in laboratory 
air and consequently in the air in balance cases is in general 
variable. A variation of O.OOOl in (O^ mole fraction is 
equivalent to a relative variation of 4 x 10 5 in calculated 
air density. Consequently, for optimum utilization of the air 
density calculation, the C() 2 abundance should be known. 
Equation (18) enables adjustment of M Q lor- departures from 
the reference level, 0.00033. of C() 2 abundance. 

6. Overall Uncertainty in p Due to 

Uncertainties in the Parameters and Variables 

in the Air Density Equation 

The random and Systematic components <>f the relative 
uncertainties in p in tables 3 and 4 when combined provide 
estimates of the components of the overall relative uncer- 
tainty in p. The random component, by quadrature, is ±3 
X 10~ 4 ; the systematic component, by addition, is ± I X 
IO -4 . The relative uncertaintv due to the variation of C0 2 
abundance (4 X 10~ 5 per 0.0001 in C0 2 mole fraction) is 
necessarily not included in this estimate. 

At 293. 15K, 101325 Pa (I atmosphere), and 50 percent 
relative humidity, the estimated overall relative uncertainties 
in p correspond to uncertainties in mass in the transfer 
between platinum-iridium and stainless steel kilogram arti- 
facts (volume difference of ~8() cm) of approximately 30 //,g 
random and 10 pug systematic 

7. Air Density at Standard Conditions, p 

The air density, p , at standard conditions (P , M Qq , R, 
T , Z ), for dry air is expressed by an equation of the form 
of (11) as 



Po 



PoM ao 
RT Z 



(33) 
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By dividing (11) by (33) we arrive at 
P\ I 7 'o\ / Z \ ( M a ' 



Po 



P», 



M a , 



(34) 



l + (e- 1) 



The gas constant, /?, has been eliminated but M a has not 
unless the apparent molecular weight of the dry air in the 
standard state, M a is equal to M a . For example, if the p is 
determined by experiment, M a is not necessarily equal to 
M a . 

If an experimental value of p of sufficient accuracy were 
available in the determination of which R was not used, R 
and its associated uncertainty would be eliminated. In the 
absence of such an experimental value, as is the present 
case, the use of an equation combining (11) with (18) is 
preferred to an equation of the form of (34). 

If standard conditions are taken to be T = 273. 15K, P 
= 101325 Pa, M a() = 28.963, and relative humidity = 0, 
for which Z = 0.99940, p is calculated by (33) to be 
1.2930 kg m~ 3 . The same result is obtained iov M a = 28.964 
(i.e., for a C0 2 mole fraction of 0.00043). 

8. Air Density Equation 

M w 
By combining (11) and (31) and substituting for e, the 

air density equation developed in this work becomes 



PMa 

RTZ 



MJ 100 P 



(35) 



By substituting the Cohen and Taylor value of/?, 8314.41 
(J A" 1 kmol" 1 ), the value 18.0152 for M w and the value 
28.963 for M^in (18), (35) becomes 



p = 0.000120273 



PMa 

TZ 



1- 1 



18.0152\ U fe s 



M„ 



100 P 



(36) 



wher 



M a 



28.963 + 12.011 [x co - 0.00033]. (37) 



For T = 293. 15K, P = 101325 Pa, 50 percent relative 
humidity and M a = 28.964 g mol -1 , the air density calcu- 
lated using (36) is 1. 1993 kg m -3 . 



8.1 . Use of Constants in the Air Density Equation 

In this section we shall investigate the effect on mass 
comparisons of the use of appropriate constant values of/, Z 
andM a in (36). 

The buoyant effect of the displacement of air by a mass 
artifact is proportional to the density of the air, p, and the 
displacement volume, V m , of the artifact. We define here the 
buoyancy correction, m b , to be added to the observed mass, 
by the following equation: 



m b = pV m = p 



(38) 



where m and p m are the mass and density, respectively, of 
the artifact. 

The variation, A/n b , in m b due to a relative uncertainty, 
Ap/p, in air density can be written 



m / Ap\ 
b.m b = — — p. 

Pm\ P) 



(39) 



By substituting a nominal value of p, 1.2 kg m 3 , in (39) 
we arrive at 



Am 6 = 1.2 



Ap 



(40) 



We shall return to this equation to estimate the uncertainty 
in mass comparisons due to the use of certain constants in 
(36) for the calculation of air density. 

To estimate the effect of a variation in / about a selected 
constant value, nominal values of P, T, M a , Z, relative 
humidity and e s are substituted in (36) and the resulting 
equation is differentiated with respect to /. Nominal values 
appropriate to the Mass Laboratory of NBS are P — 100000 
Pa, T = 298. 15K, M a = 28.964 g mol" 1 , 30 percent relative 
humidity and e s = 3169 Pa. The resulting equation is 



p = 1.1688 (1 - 0.003594/) 



Ap 
P 



0.0036 



/ ' 



(41) 



(42) 



The values of/ range from 1.0030 to 1.0046, therefore the 
maximum A// for a nominal value of/ of 1.0042 is equal to 
1.2 X 10" 3 . The corresponding Ap/p is equal to 4.3 X 10" 6 . 

To estimate the effect of a variation of Z about a selected 
constant value, from (36), 
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Ap _ -AZ 
P ' Z ' 



(43) 



To estimate the expected variation in Z, it is necessary to 
estimate the ranges of the environmental variables. A\ T and 
relative humidity. The maximum and minimum of the atmos- 
pheric pressure at NBS, taken from climatological records of 
the National Weather Service of the National Oceanic and 
Atmospheric Administration and adjusted to the elevation of 
the Mass Laboratory, are L03850 and 96160 Pa. The 
expected ranges in temperature and relative humidity are 
taken to be 18 to 28 °C and 10 to 50 percent, respectively. 
For the nominal values given in this section, the expected 
variations in the environmental variables result in a variation 
of Z between 0.99957 and 0.99971. These extremes in Z 
correspond to relative variations in Z about the nominal 
value, 0.99960 of -9 X 10 5 an d 5 X 10 \ respectively. 

The corresponding — s are 9 X 10 ' and — .-> X 10 '. 
P 

As was shown in section H.2, a variation of 0.0001 in (X) 2 
mole fraction corresponds to a relative variation in M a of 4 
X 10 5 and to a Ap/pof 4 X 10 5 . 

Combining by quadrature the maximum values of Ap/p 
due to expected variations in / and / and variation of 0.0001 
in C0 2 mole fraction, the result is 9.9 X |0 5 . This is the 
value to be used in CM)). Since the values combined arc 
maximum values, the combined value can be considered 
loosely to approximate 3 standard deviations. 

For various materials of interest in the comparison of mass 
artifacts, A////, in micrograms (/xg) has been calculated usin<: 
(39) and I kg and 100 g artifacts and tabulated in table 5. 

Table 5. Variation in liu- buoyancy correction, A///,,, resulting from (he use 
of constant values off / and l/ a in (36) 



Materia] 


Am 6 (/Ag) 


Am ft (/Ag) 




1kg 


100 g 


Platinum-indium 


5.5 


0.55 


Stainless steel 


15 


1.5 


Silicon 


54 


5.4 


Brass 


14 


1.4 


Aluminum 


44 


4.4 


Tantalum 


7.1 


0.71 


Water 


119 


12 



The variation in the apparent mass difference between 
mass artifacts due to the use of constant values of/, / and 
M a is equal to the difference in A/tif, for the artifacts. For 
example, for a comparison of kilogram artifacts of platinum- 
iridium and stainless steel the variation in the apparent 
difference in their masses is 15 — 5.5 = 9.5 fig. 

On the basis of the values calculated using table 5 and the 
precision of the balance used for mass comparisons, a 



judgment can be made concerning the adequacy of the use of 
constant values of/, Z and, M a in (36). The precision of the 
balance used for the comparisons of 1 kg mass artifacts at 
NBS is 25 /xg at the 1 standard deviation level. Thus it can 
be concluded, for example, that constant values can be used 
in the comparison of platinum-iridium and stainless steel 
artifacts, and stainless steed and brass kilogram artifacts, but 
not for comparison of platinum-iridium and silicon kilogram 
artifacts. The precision of the kilogram balance at the Bureau 
International des Poids el Mesures (BIPM) is 1.5 /Ltg at the 1 
standard deviation level, therefore it could he concluded that 
the use of constants would not be appropriate in mass 
comparisons made using that balance. 

For mass comparisons in the Mass Laboratory of NBS for 
which the values in table 5 indicate that the use of constant 
values of/ < 1.0042), / (0.99966) and l/„ (28. 904) in (36) is 
adequate, the resulting equation is 



0.0034848 
P = =--C 



0.0037960 I e 8 ) 



(43) 



For pressure in millimeters of mercury and temperature. /. 
in °C, (43) becomes 

0.46460 , 
p = - ~(P - 0.0037000 Ue 8 ). (44) 

(t + 273.15) 

For saturation water vapor pressure, e 8 , in millimeters of 
mercury (32) becomes 

e 8 = 1. 3 140 x LO 9 exp (-5315.56/0 + 27.3. 15)). (45) 

9. Conclusions 

The following are recommendations concerning the trans- 
fer of the mass unit at the various national standards 
laboratories: 1) eq (36) should be adopted for use for all 
national standards laboratories to provide both uniformity 
and the best available calculation of air density; 2) the CO2 
concentration should he treated as a variable and at least a 
"background" value should be determined for each of the 
laboratories; 3) instrumentation and practices representing 
the state-of-the-art in the measurement of the environmental 
variables should be applied. 
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